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Abstract. The sensitivity of imaging atmospheric Cerenkov telescopes (lACTs) in 

TeV 7-ray observations reachs its maximum at small zenith angles {Q < 30°) which 

provide the minimum attainable energy threshold of an instrument. However, for a 

^ I specific telescope site a number of 7-ray sources, or source candidates, can only be 

observed at much larger zenith angles (0 < 60°). Moreover the observations at large 
^^ , zenith angles allow to extend the observation time window for any object seen at small 

^^ I zenith angles, as well as to enlarge the dynamic energy range of an instrument towards 

^i^ ' the highest observable energies of 7-rays. Based on Monte Carlo simulations we present 

Qv. I here the results on the sensitivity of a stereoscopic system of 5 lACTs in observations 

"^ ' at large zenith angles. We point out some important parameters of the telescope design 

O , . which could substantially improve the efficiency of such observations with forthcoming 

' ; lACT arrays hke CANGAROO HI, HESS and VERITAS. 
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1. Introduction 

Ground-based very high energy (VHE) 7-ray astronomy explores the energy range from 
200 GeV up to 20 TeV. Recent exciting detections and observations of a number of 
7-ray sources have demonstrated the very high sensitivity of the imaging atmospheric 
Cerenkov technique (Weekes et al. 1997; Ong, 1998). The performance of this technique 
appears to be strongly dependent on the zenith angle range covered while tracking an 
object. The accessible zenith angle range is simply determined by the latitude of the 
detector and the celestial coordinates of a particular object. The 7-ray observations are 
likely to be made up to zenith angles 30° in order to detect high quality two-dimensional 
angular images of Cerenkov light from air showers. Extension of observations to larger 
zenith angles (up to 60°) might substantially widen the observational time window for 
a number of objects. Apparently that is very important for multiwavelength campaigns 
involving ground-based and satellite-born instruments in simultaneous observations of 
variable 7-ray sources, e.g. BL Lac objects. In addition observations at large zenith 
angles favour the detection of the high energy 7-rays, E > 10 TeV, which, at present, 
can be registered only using the ground-based Cerenkov technique, and set important 
constrains on 7-ray emission mechanisms. 

Sommers & Elbert (1987) first noticed that the effective detection area for 7-ray 
air showers dramatically increases in observations at large zenith angles using Cerenkov 
detectors. Expedient use of large zenith angle observations using imaging Cerenkov 
technique was suggested in the proposal of the HEGRA lACT array (Aharonian et al., 
1989). Hillas & Patterson (1991) made first simulations of Cerenkov light images from 
air showers with an inclination of 60°. The Whipple group has developed the analysis 
technique for a single stand-alone 10 m telescope (see e.g. Krennrich et al. 1999) and 
tested it by observations of Mkn 421 and Mkn 501. Large zenith angle observations 
have been used with the 3.8 m Cerenkov telescope operated by the CANGAROO group 
(Tanimori et al., 1998) to measure the energy spectrum of the Crab Nebula up to 
~ 50 TeV. Here we present results from Monte Carlo simulations of inclined air showers 
(zenith angle of 60°) for a stereoscopic system of five imaging air Cerenkov telescopes 
(lACTs). We discuss the major change in topology of Cerenkov light emission from air 
showers at large inclinations, which determines the detection rates as well as the ability 
to classify images of 7-rays. We have compared the sensitivity of the lACT system at 
small (20°) (SZA) and large (60°) (LZA) zenith angles. Current observations with the 
HEGRA lACT system generally confirm the Monte Carlo predictions (these results will 
be published elsewhere). 

2. Air shower simulations 

The ALTAI Monte Carlo code (Konopelko et al. 1996; 1999) was used to generate 7-ray 
and cosmic ray induced air showers. This code includes a detailed model of Rayleigh 
(molecular) scattering, aerosol (Mie) scattering and ozone absorption of Cerenkov 
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light in the atmosphere (Elterman, 1968; Driscoll, Vaughan, 1978). This model has 
been checked against the available experimental data for La Palma, Canary Islands 
(Hemberger, 1998). The calculations have been done with allowance for the effect of the 
geomagnetic field. By additional test simulations we found that the effect of multiple- 
scattering of Cerenkov light photons in the atmosphere is negligible because of a low 
aerosol content above the observation level (2.2 km above sea level) as well as due to 
the specific shape of the aerosol scattering-phase function (Driscoll, Vaughan, 1978) for 
SZA as well as for LZA. 

Simulations were performed for the setup of five lACTs of the HEGRA stereoscopic 
system (Aharonian et al. 1999). Each of five telescopes consists of a 8.5 m^ reflector 
focussing onto a photomultiplier tube camera. The number of photomultipliers in 
the camera was 271, arranged in a hexagonal matrix covering a field of view with a 
diameter of 4.3°. A telescope camera was triggered when the signal in two of the 
271 photomultiplier tubes exceeded a threshold of 8 photoelectrons, and the event was 
saved when at least two telescopes were triggered by Cerenkov light from air shower. 
The overall efficiency of the photon-to-photoelectron conversion was about 0.1. 

3. Lateral distribution of Cerenkov light photons 

The longitudinal profile of the number of secondary electrons in air showers of a fixed 
primary energy, E, remains almost the same at different shower inclinations when the 
distance along the shower axis is measured in gr/cm^. The depth of shower maximum 
can be calculated as T^ax = Xq • ln(E/Ec) (Xq and Ec are the radiation length and 
critical energy in air, respectively) (e.g., Nishimura, 1967). Thus a 3 TeV 7-ray 
air shower at zenith has, on average, a maximum in electron number at a depth of 
Tinax — 400gr/cm^, which corresponds to Hmax — 7.5 km above the sea. For large 
inclinations the total atmospheric optical depth increases substantially and the shower 
develops entirely in the upper layers of the atmosphere. For air showers with large 
inclinations, the height of shower maximum is far above the observation level, and the 
geometrical distance from the shower maximum to the observer is correspondingly very 
large. Thus a 7-ray-induced air shower of 3 TeV with 60° inclination on average has an 
electron maximum at Hmax — 12 km above the observation level. The corresponding 
geometrical distance from the shower maximum to observer, Lmaxj for such a shower 
is about 24 km. These geometrical factors determine the distribution of Cerenkov 
light in the plane perpendicular to the shower core. Air showers at the zenith give 
a high Cerenkov photon density in the region close to the shower core (R < 100 m) (see 
Figure 1) because of high photon emission from the low energy electrons (Ec < 1 GeV) 
which deeply penetrate into the atmosphere and suffer multiple Coulomb scattering. In 
air showers with large inclinations these electrons contribute over a much wider range of 
distances from the shower core because of the larger geometrical distance to the shower 
maximum, and as a result the corresponding mean photon density is substantially lower. 
In addition, the absorption of Cerenkov light photons in the atmosphere increases for 
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Table 1. Detection rates of 7-ray and cosmic ray- induced air showers for inclination 
of 20° and 60° and corresponding signal to noise ratio for all triggered events and for 
the high energy events above 20 TeV. 





Zenith angle: 


20° 


60° 




R7, [hr-i] 


100 


43 


All events 


RcR, [hr-1] 


5.4-104 


2.7 -lO'^ 




S/N, cr 


0.3 


0.18 




R^, [hr-i] 


1.2 


5 


E > 20 TcV 


RcR, [hr-1] 


8-102 


1.6-103 




S/N, cr 


0.03 


0.09 



air showers with large inchnations due to the large optical depth on the way from the 
point of their emission to the detector. 

The characteristic hump in the lateral distribution of Cerenkov light photons, 
caused by the emission of energetic electrons (Eg > 1 GeV) around the shower maximum, 
is shifted to larger distances from the shower axis. One can estimate the position of the 
hump using the expression Rq ~ L^ax ■ tg6'c, where 9c = 6'c(Hjnax) is the Cerenkov light 
emission angle at the corresponding shower maximum height, Hmax- Thus for a 3 TeV 
7-ray shower at the zenith and for 60° inclination the offset of the hump is 90 and 240 m, 
respectively. These "toy model" estimates are in good agreement with the simulations 
(see Figure 1). The hump becomes even more prominent at 60° inclination because of 
the reduced density of Cerenkov light photons emitted by low energy electrons. 

4. Collection areas and detection rates 

Despite the sophisticated trigger logics, the trigger condition for each telescope of a 
system ultimately relies on the size of the image (i.e., total number of ph.-e.) in the 
camera. Thus the trigger efficiency for each telescope, as well as for the entire system 
roughly reflects the lateral distribution of Cerenkov light photons at the observation 
level. Finally, the effective collection area is governed by the lateral distribution of 
Cerenkov light photons (for details see Aharonian et al., 1995). Results of calculations 
for a 5 lACT system are shown in Figure 2. A sharp increase in collection areas at 
low energies, caused by increasing trigger efficiency within the plateau of Cerenkov 
light density (R > 130 m), changes above S^ ~ 5 ■ 10*^ m^ to a logarithmic growth at 
the exponential tail of the Cerenkov light lateral distribution (see Figure 1). At large 
inclinations low energy 7-rays (E < 1 TeV) cannot trigger the telecope system because 
of the very low average image size. At the same time the broad lateral distribution 
for high energy 7-rays (E > 3 TeV) provides large collection areas which could even 
substantially exceed the collection areas at small zenith angles (see Figure 2). Thus, for 
10 TeV the 7-ray collection area for LZA is larger by a factor of 3.5 than for SZA. 

Assuming an energy spectrum of 7-rays, e.g., dJ^/dE oc E~"; a = 2.5, and a certain 
flux normalization, J7(> 1 TeV) = 10^^^ cm^^s^^, one can calculate the detection rates 
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Table 2. Effective energy tliresliold of 7-ray induced air showers at different inclination 
angles. 

e, deg 20 30 45 60 



Eth, TeV 0.5 0.7 0.9 1.8 5.0 



of 7-ray-induced air showers 



R,{>E^) = j^^S,{E)^dE (1) 



and cosmic ray showers 



dJcR 



Rcr{>E,)= dQ ScR{E,Q)-^dE (2) 

JQo Jeo dE 

above the energy threshold of Eg, where Qo is the sohd angle of the isotropic cosmic 

ray air showers, and S^ and Scr are the collection areas for 7-rays and cosmic rays, 

respectively. E is a reconstructed energy of 7-ray and cosmic ray induced air showers 

using the procedure tuned for the 7-rays (see Konopelko et al, 1999). The integral 

detection rates calculated for SZA and LZA are shown in Figure 3. Note that for 7-rays 

above 20 TeV the integral detection rate is about 4 times higher for LZA than for SZA, 

whereas the corresponding integral cosmic ray rate is higher by a factor of 2, only. Thus 

at the trigger level LZA observations reveal a substantial advantage in the detection rate 

of high energy 7-rays as well as a high signal to noise ratio, S/N = R7/(2 ■ Rcr)^''^ (see 

Table 1) ||] The signal to noise ratio for the high energy events (above 20 TeV) strongly 

depends on the energy spectrum index of 7-rays. In case of a flat energy spectrum 

{a ~ 2.0) the advantage at LZA becomes even more prominent (see below). The effective 

energy threshold of the telescope system, defined as the energy at which 7-ray detection 

rate reaches its maximum for the differential energy spectrum dN^/dE ~ E~^'^, is ~ 0.5 

TeV at zenith, and increases at larger inclinations (see Table 2). 

5. Orientation and shape of Cerenkov light images 

The above mentioned features of the shower development at large inclinations determine 
the topology of Cerenkov light images. We show in Figure 4 Cerenkov light images 
calculated for 3 TeV 7-ray-induced air shower with various inclinations. In general the 
images at large inclinations (a) contain less photons; (b) become smaller in size; (c) 
shrink to the camera center; and, (d) have a circular shape. Using our "toy model" 
considerations we may estimate the position of the image centroid (maximum of image 

I In TeV 7-ray observations the significance of 7-ray signal is calculated as S/N = ON — OFF/(ON -t- 
OFF)^/^ where ON is the number of events registered while tracking the source and OFF is the 
corresponding number of background events, Ncr. ON - OFF gives the number of detected 7- 
ray showers, N-y. In the case when OFF >> ON — OFF one can modify the formulae as follows: 
N^/(2 • Ncr)^/^ = 1^7/(2 • Rcr)^^^ • Vt, where t is the observational time and R^, Rcr are the rates of 
7-rays and cosmic rays, respectively. 
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Table 3. Acceptances of 7-ray showers after angular cut for two inclinations, 20° 
and 60°. 9 is the angular distance from the reconstructed to the actual 7-ray source 
position. 



Zenith angle: 20° 60° 

< 0.3° 0.90 0.52 

9 < 0.55° 0.95 0.70 



Table 4. Acceptances of 7-ray and cosmic ray showers and corresponding enhacement 
(Q-factor) after applying mean scaled Width cut. Air showers were simulated at two 
inclinations, 9 = 20° and 60°. 







Cut: 


«7 


thcr 


Q-factor 


9 = 


= 20° 


<w>< 1.0 
<w>< 1.3 


0.53 
0.99 


0.01 

0.15 


5.3 
2.6 


9 = 


= 60° 


<w>< 1.0 
<w>< 1.3 


0.45 
0.90 


0.05 
0.30 


2.0 

1.6 



intensity set by the emission from shower maximum) in a telescope focal plane as 
60 ~ l/tg(Ro/Linax)- At 100 m impact distance from the shower axis the centroid 
position is at ~ 1° and ^ 0.23°, for air showers at the zenith and for 60° inclination, 
respectively. The small angular size of Cerenkov light images at large inclinations 
was noticed by Hillas & Patterson (1990). Air showers with 60° inclination are very 
far from the observer which is why the image angular size is getting small in both 
longitudinal and transverse directions. The lateral spread of electrons (po) in the 
maximum of multi-TeV 7-ray shower is ~ 20 m (see e.g., Hillas, 1996). Thus the 
angular size of the Cerenkov light image along the minor image axis can be estimated 
as w ~ l/tg(Ro/Linax) — l/tg((Ro — Po)/Lmax)- K is of ~ 0.2° for shower at zenith 
and ~ 0.06° for the shower with 60° inclination (see Figure 4). To measure the 
orientation and shape of these images one needs a relatively small pixel size (angular size 
of PMTs in camera) ~ 0.1 ^ 0.15°. Note also that light smearing by optical errors may 
significantly distort the angular shape of these images. The ratio of standard second 
moment parameters. Width/ Length, shifts to larger values (see Figure 6) for showers 
at large inclinations. It shows that these images have a circular shape rather than an 
elongated elliptic shape. Consequently the image orientation is determined with larger 
uncertainties. 

For the image orientation one can use the standard A/p/ia-parameter which defines 
the angle between the major axis of an elliptical image and the line connecting the 
position of the image maximum to the camera center (we assume that telescopes are 
looking directly onto the source of 7-rays). The distribution of A/p/ia-parameter for all 
triggered telescopes in a system observing the 7-rays at 20° and 60° inclination is shown 
in Figure 5. The distribution at large inclination is relatively broader as expected. This 
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implies that 7-ray acceptance after applying the angular cut is less at large inclinations 
(see Table 3) and the corresponding enhancement factor is lower. Note that an imaging 
camera with small pixels (~ 0.1-7-0.15°) may substantially improve the angular resolution 
at large inclinations. 

In order to utilize simultaneously several Cerenkov light images for an individual 
shower the so-called mean scaled Width pararaeter (Konopelko, 1995; Daum et al., 1996) 
can be effectively used for a system of lACTs. To compensate the dependence of the 
image shape on primary shower energy and distance from shower core to the telescope 
(impact parameter) the standard parameter Width (see Fegan, 1997) (w^), calculated for 
each telescope, is scaled according to the Monte Carlo predicted values, < w >'^j, taken 
for the corresponding bin of reconstructed distance from the telescope to the shower 
core (i) and for the corresponding bin of image size (total number of photoelectrons in 
the image) (j). The mean scaled Width parameter is defined for each individual shower 
as follows 

N 

<w>= 1/N Y^w^l <w >J,. (3) 

fc=i 

where A^ is the number of triggered telescopes. The optimum cut on mean scaled 

Width is about 1.0, which gives a 7-ray acceptance of ~ 50%. However, for a precise 

determination of 7-ray spectra, a loose cut on mean scaled width (< w >< 1.2) has 

been so far used in data analysis (Aharonian et al, 1999) in order to maximize the 7- 

ray acceptance and to minimize systematic error related to cut efficiencies. In Table 4 

we show the acceptances and efficiencies of a 7-ray classification using the mean scaled 

Width cut. The small angular size of both 7-ray and cosmic ray induced air showers 

prevents effective rejection using < w >-parameter for LZA. The resulting enhancement 

(Q-f actor) does not exceed 2.0 whereas at small zenith angle it is more than 5.0. We 

may conclude that the standard orientation and shape cuts, usually used for the lACT 

system analysis, allow less significant rejection of the cosmic rays. In order to improve 

the cosmic ray rejection we introduced an additional parameter, mean scaled Length, 

< I >, defined by analogy with < w >. Two parameters, < w > and < I >, can be used 
for calculating a Mahalanobis distance, MD (see Mahalanobis, 1963), in 2-dimensional 
space 

MD = ((1- < w >)V4^> + (1- < 1 >)V<i>)'/' (4) 

where 0"<ts> and cr^j^ are standard deviations for the corresponding distributions of 

< w > and < / >. We found that the optimum value for the MD cut for LZA is 1.5. This 
analysis improves the enhancement factor by ^ 30%. Note that multivariate analysis 
technique for a single Cerenkov telescope has been discussed before by Aharonian et al. 
(1991); Hillas & West (1991). In addition we have applied the standard algorithms of 
the impact distance and energy reconstruction (e.g., Konopelko et al., 1999) used for 
the system of 5 lACTs to events simulated at large inclinations. A summary of a system 
performance at LZA is presented in Table 5. 
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6. Sensitivity to 7-ray fluxes 

The sensitivity of tlie instrument can be cliaracterized using tlie so-called sigma-per- 
hour parameter. It corresponds to the signal to noise ratio which one can achieve 
within one hour observations of 7-ray source assuming a certain DC 7-ray flux, e.g., 
Jy(> 1 TeV) = 10^"*^^ [photons/cm^s] and the energy spectrum dJ.y/E oc E^"^'^. Thus 
the sensitivity of the 5 lACT system is about 6(x and ~ la (loose cuts) at the small (20°) 
and large (60°) zenith angles. These sensitivities can be explained by the difference in 
the energy threshold and 7-ray rate for LZA and SZA. Thus, for a proper comparison of 
system performance one can use a "sigma-per-hour" for energy bin, a/hr/log{E). This 
parameter corresponds to the signal to noise ratio which one could expect detecting 7- 
rays in a fixed energy range. In present analysis we choose the corresponding energy bin 
as AEi = Ei+i — Ei, Ei+i/Ej = 1.38 which roughly corresponds to the energy resolution 
at LZA. The number of 7-ray and cosmic ray induced showers detected within the energy 
bin AEi can be calculated as 

and the corresponding signal to noise ratio is 

Si = N;^/(2-Np^)^/2_ (g) 

where (dR'^/dE) and (dR'-'^/dE) are the differential detection rates of 7-ray and cosmic 
ray air showers. 

The cr/hr/log{E) as a function of primary 7-ray energy, S = S(E), is shown in 
Figure 7. Our calculations demonstrate that for a 7-ray source with a spectrum index 
of 2.5, LZA observations have an advantage at trigger level before imaging analysis. 
After applying the software analysis cuts, the sensitivity is almost the same at 10 TeV 
with a slight advantage for LZA for higher energies. For flat-spectrum 7-ray sources 
{a ~ 2.0) the siganal-to-noise ratio is a factor of 3 higher at ~ 20 TeV for LZA compared 
with SZA. For the data shown in Figure 7 we have applied the loose analysis cuts. 
These cuts keep most of the 7-rays (/t^ ^ 0.9), which is preferable for the spectrum 
studies. However these cuts give only a modest cosmic ray background rejection with 
the corresponding quality factor of ~ 12 for SZA and ~ 5 for the LZA. Note that in 
order to achive the maximum signal-to-noise ratio one can use the tight cuts {6 < 0.22°, 
< w >< 1.0) which provide a 7-ray acceptance of /t^ ~ 0.4 and corresponding Q-factors 
as high as ^ 35 and ~ 10 for SZA and LZA, respectively. 

7. Conclusions 

We have studied the efficiency of large zenith angle observations using the stereoscopic 
system of 5 imaging atmospheric Cerenkov telescopes by means of detailed Monte Carlo 
simulations. LZA observations give a very large effective collection area at large zenith 
angles but a modest ability for 7-ray classification. The results of simulations show that 
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Table 5. Accuracy of determination of impact distance (AR), angular resolution (SQ), 
and energy resolution (AE/E), for a system of 5 lACTs. 

Zenith angle: 20° 60° 



AR, m < 10 < 50 

,56, degree 0.1 0.3 
AE/E, % 20 < 30 



for a 7-ray source with a relatively steep energy spectrum {a > 2.5) LZA observations 
provide almost the same sensitivity as normal SZA observations at the energy of 7-rays 
of about 10 TeV with some advantage at higher energies. For 7-ray sources with a flat 
energy spectrum (a^ ~ 2.0) LZA observations have a signiflcantly higher sensitivity at 
energies > 10 TeV. 

Observations at large zenith angles with the forthcoming arrays of imaging 
atmospheric Cerenkov telescopes, like CANGAROO III, HESS and VERITAS, with 
an energy threshold of 50-100 GeV, could provide an extension of dynamic energy 
range up to 20-50 TeV. Such observations would need a camera with small pixel size of 
~ 0.1 ^ 0.15° and flne optics of the telescope reflector. 
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Figure 1. Lateral distribution of Cerenkov light density in 3 TeV 7-ray- induced 
air shower with inclination angle of 0°(1); 30°(2); 45°(3) and 60°(4). The 
observation level is about 2.2 km above the sea. 
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Figure 2. Detection area as a function of a primary shower energy for 7-ray- 
induced air showers with inclination angle of 20°(1) and 60°(2). 
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Figure 3. Integral detection rates of 7-ray (a) and cosmic ray (b) air showers 
above the energy Eg for two inclination angles, 20°(1) and 60°(2). 
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Figure 4. Cerenlcov liglit images of 3 TeV 7-ray air sliower simulated at different 
inclination angles, 0°, 30°, 45° and 60° (from upper to bottom panel). The 
calculated images reproduce the two-dimensional angular distributions {Ox,Oy) 
of Cerenkov light photons hitting the telescope reflector. Angle 6x is measured 
respecting the axis which connects the telescope mirror and the shower core. 
The impact distance from the shower core to the telescope is about 100 m. 
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Figure 5. Distribution of image orientation parameter Alpha for 7-ray- induced 
air showers simulated at inclination angles of 20°(1) and 60°(2). 
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Figure 6. Distribution of Width/ Length ratio for 7-ray-induced air showers 
simulated for the inclination angle of 20°(1) and 60°(2). 
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Figure 7. Expected signal-to-noise ratio (see eqn. 5) in observations at 
small zenith angles (20°) before (curve 1) and after (curve 4) lose analysis 
cuts, assuming the differential energy spectrum dJ^/dE ~ E'^ 
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